o 
o 

> 

o 



Entangled Husimi distribution and Complex Wavelet 

transformation* 

1 2 
Li-yun Hu and Hong-yi Fan 

1 College of Physics and Communication Electronics, Jiangxi Normal University, Nanchang 330022, China 

"Department of Physics, Shanghai Jiao Tong University, Shanghai 200030, China 

November 1, 2009 



> 
in 

(N 



Abstract 

Based on the proceding Letter [Int. J. Theor. Phys. 48, 1539 (2009)], we expand the relation 
between wavelet transformation and Husimi distribution function to the entangled case. We find 
that the optical complex wavelet transformation can be used to study the entangled Husimi dis- 
tribution function in phase space theory of quantum optics. We prove that the entangled Husimi 
distribution function of a two-mode quantum state is just the modulus square of the complex 
wavelet transform of e - '" 7 ^ 2 with ip (77) being the mother wavelet up to a Gaussian function. 
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1 Introduction 



o 



Studying distribution functions of density operator p in phase space has been a major topic in quan- 
tum statistical physics. Phase space technique has proved very effective in various branches of 
physics. Among various phase space distributions the Wigner function F w (q,p) flU El C3 EJ is the 
most popularly used. But the Wigner distribution function itself is not a probability distribution due 
to being both positive and negative. To overcome this inconvenience, the Husimi distribution func- 
tion Fh (q',p') is introduced [5], which is denned in a manner that guarantees it to be nonnegative. 
On the other hand, the optical wavelet transformations have been developed which can overcome 
some shortcomings of the classical Fourier analysis and therefore has been widely used in Fourier 
optics and information science since 1980s [|6l [TJ |8l [9) . In the previous Letter [10], we have em- 
ployed the optical wavelet transformation to study the Husimi distribution function for single-mode 
case, and proved that the Husimi distribution function of a quantum state \tp) is just the modulus 
square of the wavelet transform of e~ x I 2 with if; (x) being the mother wavelet up to a Gaussian 
function, i.e., 



(lP\A h (q,p, k)\iI>) 



dxip* 



-x 2 /2 



where s = (nq + ip) , \i — y^and {tp\ Ah (q,p) is the Husimi distribution function, 



A h (q,p, k) |V>) = 2 / dq'dp'F w {q',p') exp 



, , ,2 (p' - pY 
-K{q -q) 



as well as (q,p, n) is the Husimi operator, 

A h (q,p : k) = — — : exp • 



-K(q-Q) 2 (P-Pf 



l + K 



(1) 



(2) 



(3) 
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here : : denotes normal ordering; Q = (a + a))/\/2 and P = (a — a*)/ (V2i) are the coordinate and 
the momentum operator, and ai,a\ the Bose annihilation and creation operators, [a,a^] = l,a |0) = 
0. Thus a convenient approach for calculating various Husimi distribution functions of miscellaneous 
quantum states is presented. 

Recalling that in Ref. HllH . Fan and Guo have introduced the entangled Husimi operator Ah (a, 7, k) 
which is endowed with definite physical meaning, and find that there corresponds a special two- 
mode squeezed coherent state |er, 7} representation such that Ah (a, 7, k) = |<7, 7, k) (cr, 7, «|. The 
entangled Husimi operator A^ (a, 7, k) and the entangled Husimi distribution Fh (cr, 7, k) of quan- 
tum state are given by 



and 

Fh (cr, 7, k) 



A /l (cr, 7 ,K)=4 y d2 (T 'dVA„,(cr',7')ex P |- K |cr'-a| 2 -i|7'-7| 2 }, (4) 
4 yd 2 cr'd 2 7 '^K,7')ex P |- K |a'-cr| 2 -i|7'-7| 2 }, (5) 



respectively, where F w (cr',7') = (V'l A w (cr',7') with A w (cr',7') being two-mode Wigner oper- 
ator is two-mode Wigner function. Thus we are naturally led to studying the entangled Husimi 
distribution function from the viewpoint of wavelet transformation. 

In this paper, we shall expand the relation between wavelet transformation and Wigner-Husimi 
distribution function to the entangled case, that is to say, we employ the complex wavelet trans- 
formation (CWT) to investigate the entangled Husimi distribution function (EHDF) by bridging the 
relation between CWT and EHDF. We prove that the entangled Husimi distribution function of a 
two-mode quantum state is just the modulus square of the complex wavelet transform of e I T ? I ~ / 2 
with tp (77) being the mother wavelet up to a Gaussian function. Thus we present a convenient ap- 
proach for calculating various entangled Husimi distribution functions of miscellaneous two-mode 
quantum states. 



2 Complex wavelet transform and its quantum mechanical ver- 
sion 

In Ref. H12L Fan and Lu have proposed the complex wavelet transform (CWT), i.e., the CWT of a 
signal function g (?/) by ip is defined by 

w^ g ( M , Z ) = l(iig („) r (n-i) , (6) 

f V M / 

whose admissibility condition for mother wavelets, J ^rip (rj) = 0, is examined in the entangled 
state representations (77 1 and a family of new mother wavelets (named the Laguerre-Gaussian 
wavelets) are found to match the CWT [12]. In fact, by introducing the bipartite entangled state 
representation (t] = 771 + ±772 1 , llT3l[T4l 

\V) = exp |-~ M 2 + V a\ - V *al + 0*4 J |00) , (7) 

which is the common eigenvector of relative coordinate Q\ — Q 2 and the total momentum Pi + P 2 , 

{Qi-Q2)\r))=V2r )1 \ri), (Pi + P 2 ) \rf) = V2rj 2 \t)) , (8) 

where Qj and Pj are the coordinate and the momentum operator, related to the Bose operators 

(aj,aj), [ai,aj] = 5y by Qj = (a 3 + a])/%/2 and Pj = (a- at)/(V2i) (j = 1,2), we can treat © 
from the quantum mechanically, 

W^g(n,z) = - ( f?L(i;\?L^.)(r 1 \g) = (i,\U 2 (n,z)\g), (9) 



where z — z\ + iz 2 € C, < fi G R, g (rj) = (rj\ g) , and tp (v) — ( r l\ V) are the wavefunction of state 
vector \g) and the mother wavelet state vector in (77 1 representation, respectively, and 



U 2 (n,z) 



TT 



n - 



(v\ > A* 



(10) 



aia 2 ^ In /i 




/"A 




(11) 






7T 


A*/ 





is the two-mode squeezing-displacing operator H151 1161 [TTH . Noticing that the two-mode squeezing 
operator has its natural expression in (r]\ representation B14H , 

S 2 (fi) = exp 

which is differerent from the product of two single-mode squeezing (dilation) operators, and the 
two-mode squeezed state is simultaneously an entangled state, thus we can put Eq. dlOD into the 
following form, 

U 2 z) = S 2 in) S3 (z) , (12) 
where S3 (z) is a two-mode displacement operator, S3 (z) \rj) = \q — z) and 



exp 



Pi -P 2 



V2 

= Di {-z/2)D 2 (z*/2) 



iz 2 - 



Q1+Q2 



V2 



It the follows the quantum mechanical version of CWT is 

W^g (m, C) = W S 2 <p) S3 (z) = (VI S 2 (m) Pi (-*/2) P 2 (*72) \g) 



(13) 



(14) 



Eq.([T4l) indicates that the 2D CWT can be put into a matrix element in the {rj\ representation of the 
two-mode displacing and the two-mode squeezing operators in Eq. til ID between the mother wavelet 
state vector and the state vector \g) to be transformed. Thus the CWT differs from the direct 
product of two 1 -dimensional wavelet transformations. 

Once the state vector (ip\ corresponding to mother wavelet is known, for any state \g) the matrix 
element (ip\ U 2 (/i, z) \g) is just the wavelet transform of g(rj) with respect to (ip\ . Therefore, various 
quantum optical field states can then be analyzed by their wavelet transforms. 



3 Relation between CWT and EHDF 

In the following we shall show that the entangled Husimi distribution function (EHDF) of a quantum 
state can be obtained by making a complex wavelet transform of the Gaussian function e~'' 7 ' / 2 , 
i.e., 



(ip\ Aft (<7,7,k) = e 
where \i = e x = i/k, z = z\ + iz 2 , and 



(15) 



cosh A „ / * m 

z\ = — — |7 -7-re(cr + (J)\ , 



z 2 = 



1 + K 

i cosh A 

1 + K 



[7 + 7* + K (<t — a* 



(16) 
(17) 



and (ip\ Ah (a, 7, k) is the Husimi distribution function as well as A/j (a, 7, k) is the Husimi oper- 
ator, 



A h (ct,7,k) = 



I/.- 



exp 



(ax + a\ - jj (a\ + a 2 - 7*) 



(1 + «) 



1 + K 



1 + K 



(18) 



here : : denotes normal ordering of operators. 
Proof of Eq.CE!]). 

When the transformed \g) = 1 00) (the two-mode vacuum state), noticing that (77 1 00) = e _ ' ,? ' 2 / 2 , 
thus we can express Eq. ((9]) as 



(19) 



To combine the CWTs with transforms of quantum states more tightly and clearly, using the tech- 
nique of integration within an ordered product (IWOP) H18M191I20U21H of operators, we can directly 
perform the integral in Eq. fllOD I122H 



1 [ d 2 rj 



U 2 (n, z) = - : exp 

M J 7T 



M 2 + 1, |2 , yz* + zr)* , t 



2^ l?7 



2/j 2 



/1 

sech A exp 



a 2 + a}a 2 + 77 ai — 7702 + aia2 — a[ai — a 2 a 2 — 



2/j 2 



z| z + QSjaJ tanh A + — ( z*a3, — zaj ) sech A 



x exp 



2(1 + M 2 ) 
^ajai + a5,a2^ In sech A 



z ax- za 2 . . 

exp 5 aia2 tanh A 

1 + /r 



(20) 



where we have set /i = e A , sech A = tanh A = and we have used the operator identity 

e s ° to =: exp [(e s — 1) a^a] : . In particular, when z = 0, it reduces to the usual normally ordered 
two-mode squeezing operator S 2 (/j). From Eq. d20D it then follows that 

U 2 (/x, |00) = sech A exp ( - {Zl (Zl + lZ2) + a\a\ tanh A 



2(1 + M 2 ) 



(zi — ^2) a\ — [z\ + ^2) a\ sech A > |00) 



(21) 



Substituting Eqs.([[6]), G3 and tanh A = cosh A = 5^= into Eq.lEB yields 



1 I _ , 1 2 o"y —~fcr 
3«M -^+V[/ 2 (/i,^^ 2 )|00) 



2y^ 



■ exp 



M +«|o'| KCT + 7 + 7* — K(T* + + + K — 1 I . 

(k+T) + TT^ + ^T^ - 2 + a i a 2^-r loo) - k, 7 >. 



1 + « l 2 

then the CWT of Eq . fl 1 9D can be further expressed as 



/i7T 



r? — z\ — iz 2 



ii 2 k>7>«- 



(22) 



(23) 



Using normally ordered form of the vacuum state projector 1 00) (00| —: e a i Ql a 2° 2 : .and the 
IWOP method as well as Eq.([22]) we have 



4 k 



(1 + k) 
kij* + 7 



2 : exp 



KCT i 



|7| 2 + k|ct| 2 kct + 7 f 7 
, I — , a ^ i , "2 

K + 1 1 + K 1 1 + K 



kct+7 7-kct k — 1/ + + \ + + 

H ax H a 2 H I ai a 2 + a!a 2 — a{a! — a 2 a 2 

1 + K 1 + K K+1V / 

(ai + 4 _ 7) (4 + a 2 - 7*) 



exp 



4k 

(i + «r 



k ^ai — a\ — cr^j (^a\ — a 2 — a" 

1 + K 



1 + K 



= A,, (ct,7,k) 



(24) 



Now we explain why A h (a, 7, k) is the entangled Husimi operator. Using the formula for converting 
an operator A into its Weyl ordering form H23H 

2— (-a, -f3\ A\a,P) ' exp{2 (a* ax - a\a + /3*a 2 - a\(i + a\ai + a\a 2 J } ' , (25) 

where the symbol ; ; denotes the Weyl ordering, |/3) is the usual coherent state, substituting Eq. d24D 
into Eq. (I25D and performing the integration by virtue of the technique of integration within a Weyl 
ordered product of operators, we obtain 



K 7> KK (o-,7l 



.!.()/." /■ d 2 ad 2 /3 # I ("i +4-7) (4 + a 2 - 7*) 

1 — ("«. : exp ^ 



1 + K 



/ [ a\ — 4 — cj — a 2 — cr* 

1 + K 



|a,/3> 



exp{2 ^a*ai — a|a + /3*a 2 — 4/3 + a\a± + a 2 a 2 ^} 

4' exp j — k (ai — 4 — &j (4 — a 2 — cr*^ ^ai + 4 — 7) (4 + «2 — 7*) | , 

(26) 



where we have used the integral formula 



|^ex P (cM 2 +£ 



z + 



~e-^,Re(0 <0. 



(27) 



This is the Weyl ordering form of \a, r y) Kli (a, 7] . Then according to Weyl quantization scheme H241 
|25H we know the Weyl ordering form of two-mode Wigner operator is given by 

A- w (cr, 7) = ' 5 (ai — 4 — crj <5 ^4 — a 2 — er*^ <5 + a\ — 7^ <5 + a 2 — 7*^ , (28) 

thus the classical corresponding function of a Weyl ordered operator is obtained by just replacing 



a 1 — 4 — * cr', til + 4 "~ *■ 7' ' Le -' 



4 exp — k ( ai — a 2 — (jj (^aj — a 2 — <t* J [ a\ 

4exp<{ -k|<t' - of - i | 7 ' - 7| 2 1 , 



4 - 7) (4 



+ o 2 - 7* 



(29) 



and in this case the Weyl rule is expressed as 

l°">7) KK (o-,7l = ^ J d 2 a'd 2 ^'5 (ax -a\-a}5 (a\ - a 2 - cr*) S (at + 4 - 7) 
x 5 (a\ + a 2 - 7*) [ exp j -At |er' - af - - \-y' - 



if > ■ (30) 



In reference to Eq.© in which the relation between the entangled Husimi function and the two- 
mode Wigner function is shown, we know that the right-hand side of Eq. (130D should be just the 
entangled Husimi operator, i.e. 

k 7) KK (<7, 7 | = 4 J d 2 a'd 2 7 'A w (a', 7') exp {-At |a' - a| 2 - ± | 7 ' - 7 | 2 } = A fc (<r, 7, «) , (31) 



thus Eq. Q15D is proved by combining Eqs.( 131D and 

Motivated by the proceding Letter H10H . we have futher expanded the relation between wavelet 
transformation and Wigner-Husimi distribution function to the entangled case. That is to say, we 
prove that the entangled Husimi distribution function of a two-mode quantum state \tp) is just the 
modulus square of the complex wavelet transform of e I ^ I ^ / 2 with ip (77) being the mother wavelet 

2 

up to a Gaussian function, i.e., {tp\ A h (a, 7, k) \tp) = e~*h\ 2 J J^L e -\v\ 2 /2^* ^ _ z j _ T h us 

we have a convenient approach for calculating various entangled Husimi distribution functions of 
miscellaneous quantum states. For more discussion about the wavelet transformation in the context 
of quantum optics, we refer to Refs. |[26l[27l . 
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Appendix 

We can check Eq. (13H by the following way. 

Using the normally ordered form of the two-mode Wigner operator HI lH 

A w (er,7) = — : exp {- (a\ - a\ - a^j (a\ - a 2 - er*) - (ai + 4 - 7) (4 + «2 - 7*)} : 7 (Al) 
we can further perform the integration in Eq.Q and see 

» / , . f d 2 o'd 2 i f 2 1 2 

A h (a, 7, k) = 4 / — — exp < -At \a - a\ - - |7 ~ 71 



x : exp 
4 k 



j- fai - 4 - CT ) (4 ^ a 2 - cr*) - (ai + a\ - 7) (a\ + a 2 - 7*) } : 

(at + 4 — 7 J (a\ + a,2 — 7*) At (a\ — a\ — a\ (a\ — a 2 - u*^j 



(1 + nf ' exp y i+K i+K 

= Eq.d23) =A h {a^, K ), (A2) 
which is the confirmation of Eq. (13 ID . 
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